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We will discuss here some of the recent results obtained from lattice simulations of 
QCD at non-zero temperature. Such calculations aim at a quantitative understanding of 
T— I the thermodynamics of strongly interacting matter in equilibrium. We concentrate on a 
^ discussion of the equation of state, the chiral transition in two-flavour QCD and hadronic 
^ properties related to chiral symmetry restoration in the vicinity of Tc. 

2 1- INTRODUCTION 

^ One of the most exciting predictions of QCD is the existence of a phase transition at 

^ some critical temperature to a new phase of strongly interacting matter - the quark gluon 
O plasma (QGP). Lattice calculations aim at a quantitative understanding of this phase 
2 transition and a determination of the equilibrium properties of both phases of strongly 
interacting matter. 

^ A well established result from lattice calculations is the strong dependence of the order 

of the phase transition on the QCD symmetries related to the colour as well as the flavour 
degrees of freedom. In the absence of quarks, ie. in pure SU{N) gauge theories, the phase 
transition is second order for N = 2 and first order for = 3. In the case of QCD with 
i-C n f light quark fiavours the transition is first order for > 3 and does seem to be contin- 
uous for nj = 2. The lattice calculations also show a strong dependence of the transition 
temperature on the number of partonic degrees of freedom. The transition temperature, 
expressed for instance in units of the square root of the string tension (TcZ-^/o"), is sub- 
stantially smaller in QCD with light quarks than in a pure gauge theory. In the case of 
QCD with two light quarks one finds Tc{nf = 2) ~ 150MeV, while in the purely gluonic 
theory, nj = 0, Tc is found to be substantially higher, Tc{nf = 0) ~ 260MeV. This may 
be understood in terms of a percolation threshold: In the absence of quarks there are no 
light hadrons, the lightest states are glueballs with a mass of (9(1 GeV). Therefore, a 
rather high temperature is needed to excite enough of these heavy glueball states in order 
to reach a critical density where hadrons start overlapping. 

A central goal of lattice calculations at finite temperature is to reach a quantitative 
understanding of the equation of state for strongly interacting matter and to learn details 



Oh 



*Talk given at XI. International Conference on Ultra-Relativistic Nucleus-Nucleus-Collisions, 
January 9-11, 1995, Monterey, USA. We gratefully acknowledge support through the Deutsche 
Forschungsgemeinschaft under contract No. Pe 340/3-2 



2 



30 



□ 20 



CO 



10 



0.10 



rriq/T = 0.1 







mar 



0.15 



0.20 

T(GeV) 



0.25 



0.30 



Figure 1. Energy density and pressure for two-flavour QCD on lattices with temporal 
extent = 4. 



about its structure in the vicinity of Tc, a regime which is not accessible to perturba- 
tive methods. The theoretical approach to a truly non-perturbative calculation of bulk 
thermodynamic quantities has improved a lot over recent years. In Figure 1 we show 
results from a recent non-perturbative calculation of the energy density and pressure in 
two-flavour QCD [1]. The rapid change of the energy density in a narrow region around 
Tc ~ 150MeV and the comparatively slow rise of the pressure is clearly seen. A more 
detailed discussion of the equation of state, in particular the approach to the high tem- 
perature limit and its relation to perturbative calculations is, however, hardly possible 
with the present results for QCD with light dynamical quarks. The flnite size effects are 
still quite large on the lattices presently accessible to simulations with dynamical quarks. 
Such effects have to be controlled before a comparisons with perturbative calculations 
in high temperature QCD will be possible. In the next section we will therefore restrict 
ourselves to a discussion of the equation of state in gauge theories, where a systematic 
analysis of the size dependence of lattice results is now possible [2,3]. In section 3 we 
will discuss results of studies of the order of the phase transition in two flavour QCD. A 
discussion of the influence of the chiral transition on hadron masses and other hadronic 
properties is presented in section 4. Finally we give our conclusions in section 5. 

2. EQUATION OF STATE 

The calculation of energy density and pressure in flnite temperature QCD has been a 
central target of lattice calculations. In particular, the direct calculation of the free energy 
density [4] combined with our knowledge about contributions of non-perturbative effects 
to the lattice QCD /3-function [5] allowed to eliminate remaining perturbative elements 
which were still present in earlier approaches. Nonetheless, up to very recently lattice 
calculations of bulk thermodynamic quantities have been restricted to lattices with a 
rather short extent in the temporal direction [2]. In fact, all existing results for the 
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equation of state of QCD with dynamical fermions, still are restricted to lattices of size 
X Nr with Nr = 4. 

In the following we are going to discuss a systematic analysis of the lattice equation of 
state (EOS) for pure gauge theories and its extrapolation to the continuum limit. We will 
discuss the importance of finite lattice effects for the high temperature limit of the EOS 
as well as in the vicinity of Tc. This may be viewed as a test case for future systematic 
studies of the EOS of QCD with dynamical fermions. 

2.1. The high temperature limit 

The standard forms of the lattice discretized gluon and fermion actions are known to 
lead to a strong cut-off (lattice spacing a) dependence of thermodynamic quantities on 
finite lattices of size X N^. At finite temperature these ultraviolet cut-off effects are 
controlled by aT = l/N^. Therefore they do not disappear when the "continuum" limit 
is taken at fixed N^^ ie. when the (a — > 0)-limit is mixed with the (T — > oo)-limit. These 
cut-off effects decrease with increasing Nj- like {aTY = Eor the gluonic part of the 

energy density of an ideal gluon gas (ess) one has, for instance, [2] 
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where N denotes the number of colours. On small lattices also the higher order corrections 
are still important. Eor instance, on a lattice with temporal extent Nj- = 4 finite size 
effects lead to a 50% distortion of the free gluon gas, ie. ea(4) is 50% larger than the 
corresponding continuum Stephan-Boltzmann constant. On a twice as large lattice this 
reduces to less than 10%, ea(8) ~ 1.07. It is well known from calculations on lattices with 
Nt = 4 that, for instance, the energy density above Tc rapidly approaches this "distorted" 
ideal gas value. However, in how far perturbative and/or non-perturbative contributions 
to the energy density are hidden under the large finite size corrections can only be clarified 
through a systematic analysis on various size lattices. 

In the pure gauge sector, ie. for SU (2) [2] and SU (3) [3] gauge theories, we have recently 
performed detailed calculations of thermodynamic quantities on lattices with temporal 
extent Nj- = 4, 6 and 8. In Eigures 2a and 2b we show results for the energy density and 
pressure in a SU (3) gauge theory. These have been obtained from calculations on lattices 
of size 16"^ X 4 and 32"^ x Nj- with Nj- = 6 and 8. As can be seen the finite cut-off effects 
visible in the calculation of the pressure and energy density are similar to those calculated 
in the ideal gas limit. We note that the finite size effects are largest for high temperature. 
This is to be expected as the finite cut-off mainly affects the high momentum modes, 
which dominate the energy density and pressure in that limit. Still on a quantitative level 
these effects are smaller than expected on the basis of the perturbative calculation for an 
ideal gas. This suggest that low momentum modes, which are not sensitive to the short 
distance cut-off effects, still provide a relatively large contribution to the energy density 
and pressure. 

We can use the leading 1/A^^ dependence of energy density and pressure for an ex- 
trapolation of the Nt = 6 and 8 results shown in Eigure 2 to the continuum limit. This 
suggests that already the Nt = 8 results themselves provide a very good approximation 
to the continuum limit. 
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Figure 2. Energy density (a) and pressure (b) in a SU{3) gauge theory on 16"^ X 4, 32"^ X 6 
and 32"^ X 8 lattices. The dashed horizontal lines indicate the corresponding result for an 
ideal gas on these lattices. The solid horizontal line shows the continuum ideal gas result. 
The vertical bar in (a) gives the size of the latent heat discontinuity at Tc [6]. 



We note that the increase in pressure is much slower compared to the that in the energy 
density. At T = 2Tc the energy density deviates from the ideal gas value only by about 
15% while the pressure is still more than 30% below the Stephan-Boltzmann value. This 
has strong consequences for the velocity of sound, which deviates in this temperature 
interval above Tc strongly from the ideal gas value, = 1/3 (Figure 3). 

The complicated infrared structure of QCD is expected to play an important role also 
for the thermodynamics at high temperature. Non-perturbative modifications of the low 
momentum spectrum are expected to be responsible for a substantial part of the deviations 
of the FOS from ideal gas behaviour, e = 3p. The leading ideal gas term is eliminated in 
the difference e — 3p. If infrared physics is mainly responsible for the deviations of this 
quantity from zero, we should expect that it is less sensitive to finite cut-off effects then 
e and p separately. This is indeed the case as can be seen from Figure 4. 

The difference e — 3p is given by the trace of the energy-momentum tensor normalized 
to its zero temperature value. It thus may be viewed as the definition of the temperature 
dependent gluon condensate [7] 

e-3p = {G')o-{G')T . (2) 
where (G'^)o denotes the gluon condensate at T = 0. Perturbative corrections start at 
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^ N{N'-l)g{T) + 0{g') 
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Figure 3. The velocity of sound in a SU{3) 
gauge theory. Results for the 32"^ X 8 lattice 
are indistinguishable from those shown for 
the 32^ X 6 lattice. 



Figure 4. The interaction measure e — 3p 
normalized to T"^ for a SU (3) gauge theory 
on three different lattices. 



The difference (e — 3p)/T'^ thus is expected to approach zero only logarithmically with 
increasing temperature. From the high temperature behaviour shown in Figure 4 we 
find that g{T)>l is needed to describe this regime with the leading order perturbative 
correction. Unfortunately, the perturbative expansion cannot be trusted for such large 
values of g{T) [8]. 

In order to express (e — 3p) in physical units and relate the numerical results to the 
zero temperature estimates for the gluon condensate we have to set the scale for the 
temperature. This can, for instance, be achieved by using results from lattice QCD for 
the ratio Tc/y^ with ~ 420MeV. In order to do so we have re-calculated the critical 
couplings for the deconfinement transition on 32"^ X Nj- lattices with Nj- = 6, 8 and 12. 
At these couplings we also have calculated the string tension on lattices of size 32"*. From 
this we find Tjy^ ~ 0.62 [3]^. 

(e — 3p)/T'^ shows a pronounced peak close to T/Tc = 1.1. Using the result for Tc/y^ 
we deduce from Figure 4 

^-^^- 2.4^ at r = i.ir., (4) 



r4 fm^ 

which should be compared with estimates for the gluon condensate at zero temperature, 
{G'^)o = 2.0 Gev/fm"^ [7]. This suggests that the constant zero-temperature gluon con- 
densate gives the dominant contribution to (e — 3p)/T'^ close to Tc, ie. {G'^)t turns out 



^This is about 10% larger than our earlier estimate [9], which was based on old results for the critical 
couplings calculated on rather small spatial lattices. 
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to be small at T = l.lTc. It becomes negative in this temperature region and approaches 
the negative perturbative behaviour indicated by Eq. 3. 

2.2. The critical energy density 

In the vicinity of the critical temperature the dependence of the energy density on Nj- 
itself is small. Standard renormalization group arguments (finite size scaling) suggest, that 
in this region the dominant size dependence is not related to the non-vanishing lattice 
spacing but rather to the large correlation length close to Tc, ie. the low momentum 
structure of the free energy. In fact, the behaviour of thermodynamic quantities in this 
region is controlled by the singular part of the free energy density, 

fsit,V) = ^^Q{tiTV'^y^^) , (5) 

where t = [T — Tc)/Tc denotes the reduced temperature, Q is an arbitrary scaling function 
and u denotes the correlation length critical exponent. In terms of the lattice parameters 
N(j and Nj- we have = NJN^. Finite size effects in the vicinity of Tc are thus 

controlled by this ratio rather than by aT = l/Nj- alone. The scaling behaviour in terms 
of TV^^^ has clearly been seen in the case of an SU (2) gauge theory and could be used 
there to extract the critical energy density [2] in the thermodynamic limit. The SU{2) 
gauge theory has a second order deconfinement phase transition. The critical behaviour 
is controlled by critical exponents of the 3-d Ising model. One thus expects the critical 
energy density to scale like 

^, = (^) + const. (TV^/^f-^y^ , (6) 

where a and u denote the critical indices of the 3d Ising model, [a — l)/^ = —1.41. As 
can be seen from Figure 5 this scaling behaviour holds very well. We note that the critical 
energy density turns out to be rather small, [cc/T^) = 0.256(23) or ec/esB — 1/8. 

A similarly detailed analysis of the finite volume effects close to Tc does not exist for 
the energy density of the SU (3) gauge theory. However, the latent heat at the phase 
transition point has been studied on fairly large lattices [6]. The experience from the 
SU{2) analysis suggest that the energy density at Tc in the plasma phase, e'^[Tc)^ provides 
a good estimate for the critical energy density. From Figure 2 we find e'^[Tc)/T^ ~ 2. 
Using again the string tension to set a physical scale we find here a critical energy density 
of about 1 GeV/fm"^. A similar value is also obtained from the energy density calculation 
for two-fiavour QCD shown in Figure 1. We deduce from this Figure €c/Tc = 6ib4. With 
a critical temperature of 150 MeV this also yields a rather small critical energy density, 

Cc < 1 GeV/fm3 . (7) 

We note, however, that this estimate (as usual) is rather sensitive to the value of Tc and 
might still change with more accurate determinations of the chiral transition temperature. 
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Figure 5. Extrapolation of the critical energy density for the SU{2) gauge theory as 
function of TV^^^ = N^/N^ to the infinite volume limit. Shown is a collection of data 
from various lattice sizes. 



3. THE CHIRAL PHASE TRANSITION 

The nature of the chiral phase transition does seem to depend strongly on the number 
of light quark fiavours, nj. Monte Carlo calculations performed so far with two light 
fiavours suggest that the transition is continuous, while it seems to be first order for 
nj > 3. The special role of two-fiavour QCD has been noticed already some time ago by 
Pisarski and Wilczek [10]. They have suggested that the dynamics of the chiral phase 
transition is controlled by an effective, 3-c/ scalar Lagrangian constructed in terms of the 
chiral order parameter field, a ~ ?/'?/'. This effective theory has a global 0(4) symmetry 
{SU{2) X SU{2)), which is spontaneously broken. It therefore is expected that in the case 
of a continuous transition the critical behaviour is governed by critical exponents of a 
3-0?, 0(4) symmetric spin model [10,11]. A remaining question in this scenario is to what 
extent the axial U{1) influences the critical behaviour in the vicinity of Tc [11,12]. An even 
more fundamental problem is whether an effective scalar theory can at all describe the 
critical behaviour of a theory, in which the scalar flelds are only constructed as fermionic 
bilinears. It recently has been argued by Kocic and Kogut that this may not be the 
case and that two-flavour QCD may show mean held behaviour in the vicinity of Tc [13]. 
On the lattice an additional problem arises: None of the fermion Lagrangians used for 
numerical calculations has the complete chiral symmetry of the continuum Lagrangian. In 
the staggered fermion formulation, for instance, one has a chiral U{1) rather than the full 
SU{2) symmetry, which only is restored in the continuum limit. This too may influence 
the critical behaviour in numerical studies. A detailed quantitative investigation of the 
critical behaviour thus is needed. 
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Figure 6. The chiral susceptibility versus 6/(/^ for three values of the quark mass on a 
8^ X 4 lattice. 
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The critical behaviour of QCD with light quarks is controlled by two dimensionless 
parameters, the reduced temperature t = {T — Tc)/Tc and the scaled quark mass, h = 
niq/T. The latter corresponds to an external magnetic field in spin models. In the case of 
QCD with dynamical quarks it is the quark mass rather than the spatial lattice volume, 
which is the most stringent limitation for large correlation lengths in the vicinity of the 
critical point. Rather than performing a finite size scaling analysis, as it has been done in 
the case of the deconfinement transition in SU{N) gauge theories, it is more appropriate 
to consider the scaling behaviour in terms of the quark mass. Finite volume effects may 
in a first approximation be ignored. 

In the vicinity of the critical point the behaviour of bulk thermodynamic quantities 
is related to thermal (yt) and magnetic (yh) critical exponents, which characterize the 
scaling behaviour of the singular part of the free energy density, 

fs{t,h) = -^lnZ, = b-'f{byH,by-h) . (8) 

Here b is an arbitrary scale factor. Various scaling relations for thermodynamic quantities 
can be derived from Fq.8 [14,15]. For instance, one finds for the chiral order parameter, 
{tptp)^ and its derivative with respect to the quark mass (the chiral susceptibility Xm) 



{i)i)){t,h) = h^l^F{z) 



(9) 



with scaling functions F and F' that only depend on a specific combination of the reduced 
temperature (t) and scaled quark mass (/i), z = th~^^^^. The critical exponents f3 and S 
are given in terms of yt and yh as (3 = [l — yh)/yt and 5 = y^j il — yh). A direct consequence 
of Fq. 9 is, for instance, that peak located at some fixed value of Zc- This defines 
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a pseudo-critical point, tc{h) = Zch^^^^^ at which the peak height of Xm increases with 
decreasing quark mass according to Eq. 9. Figure 6 shows the scahng behaviour of Xm 
for three values of the quark mass [15]. The rapid rise of Xm with decreasing quark mass 
can be used to determine the critical exponent 5. 
A quite useful observable is the chiral cumulant, 



A(z] 



hXm 

ixx) 



(10) 



From Fq. 9 it follows that A is directly a scaling function, ie. it only depends on z. 
Moreover, it is uniquely determined at 2; = 0, where it takes on the value 1/5 even at non- 
vanishing values of the quark mass. The chiral cumulant thus allows a direct calculation of 
this critical exponent. In Figure 7a we show a collection of the pseudo-critical couplings 
obtained from Monte Carlo calculations with various values of the quark mass and on 
various different lattice sizes. The curves show various fits for the quark mass dependence 
of these pseudo-critical couplings. Apparently it is difficult to distinguish 0(4) exponents 
from mean field exponents in this way. We note, however, that the location of the zero 
quark mass critical coupling is quite well determined. The cumulant A is shown in 
Figure 7b. For the two smaller quark masses this quantity has been evaluated in the 
region of the estimated zero quark mass critical coupling, (3c{mq = 0) ~ 5.243(10). In this 
interval the cumulant takes on values between 0.21 and 0.26, which is compatible with 
the 0(4) exponent 1/5 = 0.205(45). 
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Figure 7. The pseudo-critical couplings (a) of the chiral transition in two-fiavour QCD 
on lattices with temporal extent Nj- = 4 and the chiral cumulant (b). Shown are results 
from calculations with different values of the quark mass on a 8"^ X 4 lattice. 



Further scaling relations can be obtained by considering the quark mass dependence of 
the thermal susceptibility, Xt ~ d{'tp'tp) / dt^ and the specific heat. A combined analysis of 
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these allows a determination of the critical exponents yt and yh- Results from the first 
analysis of this kind are summarized in Table 1 and compared with the known exponents 
for 3-0?, 0(4) symmetric spin models [16] and mean field behaviour. 

Table 1 



Thermal (yt) and magnetic (yh) exponents from simulations of two-fiavour QCD on an 
8"^ X 4 lattice (LGT). Also given are mean field exponents (MF) and those for a 3-c/, 0(4) 
symmetric spin model. 





LGT 


0(4) 


MF 


yt 


0.69(7) 


0.446(5) 


0.5 


Vh 


0.83(3) 


0.830(6) 


0.75 



As can be seen from the Table the exponent yh extracted from lattice calculations is in 
very good agreement with the expected 0(4) model behaviour. It corresponds to a value 
of 0.205(45) for the exponent 1/^, which is significantly below the value 1/3 expected 
for a mean field critical exponent. The exponent j/f, however, comes out too large in the 
present calculation. It may be expected that the determination of yt is more sensitive to 
finite volume effects than that of yh- In the limit of vanishing quark masses it is yt that 
controls the finite volume scaling behaviour of thermodynamic quantities. 

Clearly the analysis of critical exponents on a 8"^ X 4 lattice can only be considered as a 
first exploratory step towards a detailed investigation of critical exponents in two-fiavour 
QCD. In the near future calculations on larger lattices with smaller quark masses will 
certainly lead to much better determinations of critical exponents. This should allow to 
distinguish between mean field and 0(4) critical behaviour. 

4. CHIRAL BEHAVIOUR CLOSE TO T, 

Fffective theories, deduced from QCD in order to describe low energy properties of 
hadrons, establish a close link between hadronic properties and the non-perturbative 
structure of the QCD vacuum, which is described by various non-vanishing condensates. 
This is, for instance, used in the operator product expansion (OFF) for correlation func- 
tions of hadronic currents [17]. It allows to relate hadron masses to condensates of the 
quark and gluon fields. The temperature and density dependence of the latter has been 
discussed within the context of chiral perturbation theory [7]. The applicability of the 
OFF as well as chiral perturbation theory is limited to low temperatures. In the vicinity 
of the QCD phase transition from hadronic matter to the quark gluon plasma phase, how- 
ever, the hadronic medium becomes so dense that the infiuence of many particle states 
and resonances can no longer be neglected. In this regime an entirely non-perturbative 
approach becomes mandatory. 

The numerical results for the energy density discussed in Section 2 suggest that up to 
temperatures of about 0.9Tc the hadronic medium is well described by a gas of hadrons, 
in which only the lightest states contribute. The slow variation of the energy density 
and pressure below 0.9Tc shows that also the gluon condensate will vary only little with 
temperature. This is also expected from chiral perturbation theory, which suggests that 
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the temperature dependence of {G'^)t starts only at 0{T^) [7]. For the chiral condensate 
(?/'?/') chiral perturbation theory gives, however, a stronger dependence on temperature, 

n2 - 1 / \ n2 - 1 / \^ . «J 



(T) /]/)!/) (T = 0) 



o N^=0(N^ = 8) 
V N^=0(N^ = 4) 
« N^=2(N^ = 4) 
□ Nj^=2 + 1 (N^=4) 
A Nj,=4(N^ = 8) 
0.5 



T/T 



Figure 8. The normalized chiral condensate versus T/Tc. Shown are results from sim- 
ulations with various numbers of flavours. Only in the pure gauge case (nj = 0) have 
the results been extrapolated to vanishing quark mass. The zero quark mass extrapola- 
tion for QCD with dynamical quarks may still lead to a somewhat stronger temperature 
dependence. 



In Figure 8 we show a collection of results for (tptp) obtained from lattice calculations 
for QCD with various number of flavours. They suggest at most a weak dependence of 
the chiral condensate on temperature, which may already indicate that we cannot expect 
to find drastic modifications of hadronic properties at least up to temperatures T ~ 0.9Tc. 

Information about hadron masses and decay constants can be extracted in lattice sim- 
ulations from the long-distance behaviour of correlation functions of hadron operators 

Gh{x) = {H{x)H+{0)) ^ e-''^"^^^ , x = {t,x) , (12) 

with H[x) denoting an operator with the appropriate quantum numbers of the hadronic 
state under consideration, for instance H[x) = 4'u{x)^5tpd{x) for the pion. At zero tem- 
perature one studies the behaviour of the correlation function at large Fuclidean times r. 
From the exponential decay of Gh at large values of r one deduces the mass of the lightest 
state in this channel, whereas the amplitude is related to the decay constant of such a 
hadronic state. At finite temperature the Fuclidean time extent is limited, < r < 1/T. 
It thus is preferable to study the behaviour of Gh for large spatial separations, \x\ — > oo. 
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This yields information about the finite temperature screening masses which are related 
to pole masses as long as there is a bound state in the quantum number channel under 
consideration. 

Finite temperature screening masses have been studied in lattice simulations of QCD 
for quite some time [18]. A very drastic qualitative change in the screening masses is seen 
when one crosses the QCD transition temperature. Parity partners become degenerate 
above Tc, the pseudo-scalar mass becomes massive and approaches the free quark/anti- 
quark value, mineson = 27rr, at large temperatures. These features do not seem to depend 
much on the number of quark fiavours. In particular, they also have been found in 
quenched QCD simulations. It thus seems to be meaningful to first study the behaviour 
of hadronic properties in the quenched approximation where results on large lattices with 
high accuracy can be obtained. In the following we will describe the results of such an 
investigation performed on a rather large lattice (32"^ X 8) close to the phase transition 
temperature [19]. 

4.1. The GMOR Relation and the Pion Decay Constant 

The GMOR relation relates the chiral condensate to the pion mass and pion decay 
constant, 

f^ml = mq{tptp)^^=Q . (13) 

Below Tc the pion is a Goldstone particle, its mass squared depends linearly on the quark 
mass, 



(14) 



A calculation of the chiral condensate and the pion mass at different values of the quark 
mass allows the determination of the pion slope, a^, and the zero quark mass limit of the 
condensate, {tptp)mq=o- 

The pion decay constant can be determined directly from the relevant matrix element 

V2Am2 = (0|V'„75V'dk+) • (15) 

The square of the matrix element appearing on the right hand side of Eq. 15 is propor- 
tional to the amplitude of the pion correlation function and can be determined in a Monte 
Carlo calculation. A comparison of determined this way with the ratio {tptp)m.q=o/0'TT 
thus provides a direct test of the GMOR relation at finite temperature. This is shown in 
Figure 9a at T ~ 0.92rc. Moreover, we can compare the value of calculated at finite 
temperature with corresponding zero temperature results. This is shown in Figure 9b. 
Clearly we neither have evidence for violations of the GMOR relation nor for a significant 
change of with temperature below Tc. 

The sudden change of "/jt" above Tc refiects the drastic change in the structure of the 
pseudo-scalar correlation function. The latter does no longer have a pole, which would 
correspond to a Goldstone-particle (pion). The existence of a pseudo-scalar bound state 
above Tc thus is questionable. Certainly for large temperatures such a state does not 
exist, the correlation function is dominated by a quark/anti-quark cut. The notion of 
used for the square root of the amplitude of the correlation function should thus be used 
with caution above Tr. 




Figure 9. The pion decay constant at T = 0.92rc as a function of the quark mass (a) and its 
temperature dependence (b). In Figure (a) we show a comparison of results in quenched 
QCD [nj = 0) obtained from the amphtude of the pion correlation function through 
Fq.(15) (circles). The zero quark mass extrapolation is compared with the result obtained 
from an application of the GMOR relation (star). Also shown is the zero temperature 
result at this value of the gauge coupling (square). 



4.2. Vector Meson Mass and Nucleon Mass 

The possibility of a variation of the p meson mass with temperature is particularly in- 
teresting because of its possible experimental consequences. It might lead to modifications 
of dilepton spectra, which are experimentally detectable. The temperature dependence of 
the meson masses has been discussed within the framework of the OFF. Arguments have 
been given that the meson masses are temperature independent up to 0{T^). The Monte 
Carlo calculations of the vector meson correlation function at finite temperature also show 
no significant temperature dependence of the mass even close to Tc- In Figure 10a we 
show the result of our calculation of the screening mass in the vector channel correla- 
tion function at T ~ 0.92rc and compare this with zero temperature calculations at the 
same value of the gauge coupling. There is no evidence for any temperature dependence. 
The same holds true for the nucleon mass, although the details are more subtle in this 
case. As can be seen in Figure 10b there is a clear difference between the local masses, 
m]\j{z) ~ \nGN{z) I Gn{z -\- 1), extracted on a large zero temperature lattice and those 
extracted at finite temperature on a lattice of size 32"^ X 8. However, as the nucleon is a 
fermion, there is a non-negligible contribution from the non-zero lowest Matsubara mode 
to the nucleon screening mass, 

rriN = ^Jm% + [ttTY (16) 

After removing the thermal contribution, ttT, the nucleon mass, mjv, agrees with the zero 
temperature result within statistical errors. 
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Figure 10. Temperature dependence of hadron masses: In Figure (a) we show a compari- 
son of /9-mass calculations at zero and non-zero temperature, T = 0.92rc for a fixed value 
of the gauge coupling, 6/(/^ = 6.0 and for various values of the quark mass. In Figure (b) 
estimates for the nucleon mass obtained from the nucleon correlation function at distance 
z are shown. These local masses converge to the nucleon mass in the limit z ^ oo. The 
horizontal band indicates the corresponding extrapolation at zero temperature. (See text 
for appropriate subtraction of the finite temperature Matsubara mode). 



5. CONCLUSIONS 

We have discussed three topics in finite temperature QCD: the structure of the equation 
of state, the order of the chiral phase transition and the dependence of hadron properties 
on temperature. A detailed quantitative understanding of these issues plays an important 
role in our attempts to clarify the structure of QCD and its non-perturbative infrared 
behaviour at high temperature. 

Studies of the equation of state for pure gauge theories have now reached an accu- 
racy, where systematic effects in the lattice approximations can be controlled and an 
extrapolations to the continuum limit can be performed with some confidence. We can 
expect that similarly accurate results will be available for QCD with dynamical fermions 
in the very near future. This then will provide important input for the description of the 
hydrodynamic expansion of hot hadronic matter created in heavy ion collisions. 

The numerical analysis of the chiral transition will help to understand the relevant 
massless degrees of freedom in the vicinity of Tc, the role of composite scalar fields for the 
dynamics of the phase transition and the interplay between the deconfinement and chiral 
symmetry restoration mechanisms. First steps in this direction have now been undertaken 
through quantitative studies of critical exponents for two-fiavour QCD. 

A calculation of the temperature dependence of hadron properties is of immediate 
relevance for the discussion of possible experimental signals for the QCD phase transition 
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in heavy ion experiments. Any change of hadron masses, their widths and decay constants 
should lead to experimentally observable consequences. Although drastic changes of the 
particle spectrum have been observed in lattice calculations, so far no sign of a temperature 
dependence has been found below Tc. This may be due to the current limitations for 
the study of hadron properties at finite temperature (small lattices, large quark masses, 
quenched approximation). This issue certainly will be studied in much more detail with 
more realistic parameters in the future. 
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